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Abstract The present work is based on a description for the angular mometum loss rate due to 
magnetic braking for main-sequence stars on the relationship between stellar rotation and age. In 
general, this loss rate denoted by dj/dt depends on angular velocity ii in the form dJ/dt oc ii q , 
where q is a parameter extracted from nonextensive statistical mechanics. Already, in context 
of stellar rotation, this parameter is greater than unity and it is directly related to the braking 
index. For q equal to unity, the scenario of saturation of the magnetic held is recovered, otherwise 
q indicates an unsaturated Held. This new approach have been proposed and investigated by de 
Freitas & De Medeiros for unsaturated Held stars. In present work, we propose a nonextensive 
approach for the stellar rotational evolution based on the Reiners & Mohanthy model. In this 
sense, we developed a nonextensive version of Reiners & Mohanthy torque and also compare this 
generalized version with the model proposed in de Freitas & De Medeiros based on the spin-down 
Kawaler torque for the main-sequence stars with F and G spectral types. We use the same sample 
of ~ 16000 field stars with rotational velocity usin i limited in age and mass. As a result, we 
show that the Kawaler and Reiners & Mohanthy models exhibit strong discrepancies, mainly in 
relation to the domain of validity of the entropic index q. These discrepancies are mainly due to 
sensitivity on the stellar radius. Finally, our results showed that modified Kawaler prescription is 
compatible with a wider mass range, while the Reiners & Mohanty model is restricted to masses 
less than G6 stars. 


INTRODUCTION. — It is generally accepted that 
the magnetic braking is the fundamental mechanism re¬ 
sponsible for angular momentum loss through magnetic 
stellar winds in several classes of stars, such as main- 
sequence field stars and cluster stars. This break was ini¬ 
tially suggested by Schatzman [l] in 1962. Later, Kraft 
[2] showed that the behavior of mean rotation velocity 
of low-mass-main-sequence stars below 1.5M 0 is preferen¬ 
tially due to magnetic wind. As mentioned by Kawaler 6], 
main-sequence (hereafter MS) stars with mass less than 
1.5 Mq contain their deep surface convection zones which 
support magnetic field as well as acoustically controlled 
stellar winds. 

Skumanich [3]’s pioneering work reveals that the an¬ 
gular momentum loss of G-type main-sequence stars for 
Hyades and Pleiades obeys the following relationship 
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dJ/dt oc fl 3 , where denotes angular velocity. Thus, the 
projected rotation velocity usini decreases roughly with 
f -1 / 2 , is valid only within the equatorial plane. Here, t 
is stellar age and exponent —1/2 is related to exponent 
from dJ/dt by simple relationship f -1 / 2 = t 1 /! 1-3 ), where 
dJ/dt is the angular momentum loss rate. In addition, un¬ 
saturated G stars are in agreement with this relationship. 
As quoted by Schrijver & Zwaan [2] , empirical Skumanich 
relationship, revisited by several authors, has stood the 
test of time, albeit that the equations fitted to the data 
do differ somewhat from author to author. In contrast, 
Chaboyer et al. [5] elaborated a general parametrization 
developed by Kawaler |6j and showed that stars that rotate 
substantially faster are subject to magnetic field satura¬ 
tion at a critical rotation rate ui sa t- Krishnamurthi et al. 
71 proposed the inclusion of a Rossby scaling at satura¬ 
tion velocity for stars more massive than 0.5 Mq. However, 
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Sills, Pinsonneault & Terndrup |8j argue that the Rossby 
number is inadequate for stars with mass less than 0.6 M@. 

Recently, Reiners & Mothanty 91 proposed a new brak¬ 
ing approach based on the strong dependence on the stellar 
radius which arises from definition of the surface magnetic 
strength (Bq oc Q a ). In contrast, Kawaler |6 assumes 
that rotation is related to the magnetic field strength 
( B 0 R 2 oc fl a ) reducing dependence on the radius. Accord¬ 
ing to 91, the factor of R 2 implies deep changes to the an¬ 
gular momentum loss law. Another important difference 
between the modified Kawaler and Reiners & Mothanty 
torques is in the saturation threshold u! sa t- Reiners & 
Mothanty [9] assume that w SQt does not depend on mass 
and it only affects the torque in unsaturated stars. Their 
work showed that the angular momentum loss law for low- 
mass stars obeys the relationship dJ/dt oc ST 5 . For this 
torque, the rotational velocity decreases with stellar age 
according to t -1 / 4 . 

More recently, de Freitas & De Medeiros 10 revis¬ 
ited the modified Kawaler parametrization proposed by 
Chaboyer et al. [5 in the light of the nonextensive sta¬ 
tistical mechanics 11|. This new statistical theory were 
introduced by Tsallis lT], and further developed by Cu- 
rado & Tsallis 12| and Tsallis, Mendes & Plastino 13 


with the goal of extending the regime of applicability of 
Boltzmann-Gibbs (B-G) statistical mechanics that present 
just plainly fail when applied in out from equilibrium sys¬ 
tems. de Freitas & De Medeiros 110 analyzed the rota¬ 
tional evolution of unsaturated F- and G- field stars lim¬ 
ited in age and mass in the solar neighborhood. They 
attribute to the entropic index q extracted from the Tsal- 
lis formalism as a parameter that describes the level of 
magnetic braking as called “braking index”. They also 
associated this parameter with the exponent of dynamo 
theory (a) and to magnetic field topology ( N ) through the 
relationship q = 1 + 4aN/3. As a result, they showed that 
the saturated regime can be recovered in the nonextensive 
context assuming the limit q —>■ 1. This limit is particu¬ 
larly important because it represents the thermodynamic 
equilibrium valid in Boltzmannian regime. Indeed, the 
torque in nonextensive version is given by dJ/dt oc SI 9 re¬ 
vealing that the rotational velocity of F- and G-type main 
sequence stars decreases with age according to t 1 ^ 1 ~ q \ 
The values of q obtained by de Freitas & De Medeiros 10 


qo is the index q for the cluster with age to- In general, 
this behavior indicates that stars loss the memory of past 
angular momentum history when their rotational distribu¬ 
tion becomes extensive, i.e., q = 1. However, we wonder if 
this result set prevails when invoke another parametriza¬ 
tion model of the angular momentum loss, in our specific 
case, the Reiners & Mothanty torque. 

Epstein & Pinsonneault 


16 compare the prescriptions 


for the strength of stellar wind of modified Kawaler and 
Reiners & Mohanty torque with and without the disk af¬ 
fect for stars with mass between 0.5 and I.OMq. They 
conclude that the largest differences between two prescrip¬ 
tions appear among low mass stars. As a result, they claim 
that the Reiners & Mohanty torque preserves rapid rota¬ 
tors of low mass stars. This effect is related to the Rein¬ 
ers & Mohanty torque predicts a slower spin down than 
Kawaler does for the first few Myr. Epstein & Pinson¬ 
neault 16 still claim that Kawaler parametrization is in¬ 


distinguishable to the Skumanich law beyond the Zero Age 
Main-Sequance (ZAMS). Indeed, this result is expected 
because the values adopted for a and N imply that angu¬ 
lar velocity decays with time asymptotically as ST ~ t -1//2 . 
In contrast, for the Reiners & Mohanty torque using same 
values of a = 1.5 and N = 2, the spin down rate de¬ 
creases to a weaker power law, i.e., SI ~ f -1 / 4 . In fact, 
the Kawaler (standard) framework of angular momentum 
loss rate described above is limited to stars below late-F 
stars. For stars above early-F a simple scaling of the solar 
torque by w 3 does not explain the fast rotators in young 
stars, in this case, saturation in the angular momentum 
loss is necessary. In both loss laws, rotation rate from sat¬ 
urated stars diminishes exponentially rapidly to critical 
saturation velocity. In general, according to de Freitas & 
De Medeiros 10 , as the exponent of ui is sensitive to stel¬ 


point out that it has a strong dependence on stellar mass. 

Some applications have been made to test the nonex¬ 
tensive approach proposed by de Freitas & De Medeiros 
10 , among them Silva et al. 14 and de Freitas et al. 


lar mass, we obtain fluctuations this exponent in narrower 
mass range than those above defined. 

In this study we compare the Kawaler and Reiners and 
Mothanty torques in the context of the Tsallis’ nonexten¬ 
sive statistical mechanics. In Section 2, we revisit and 
generalize the parametric Reiners and Mothanty model 
for angular momentum loss by magnetic stellar wind, as 
well as, we compares both models. Finally, conclusions 
are presented in Section 5. 

Generalized Reiners &: Mohanty torque. — Rein¬ 
ers & Mothanty | : 9 propose a new approach for the stellar 
magnetic braking scenario. This approach considers two 
major deviations from modified Kawaler prescription (see 


15 using this model for open cluster stars. In Silva et al. 
14 , the authors define the index q extracted from mod¬ 
ified Kawaler parametrization as qx , where the subscript 
K is an abbreviation for Kawaler. In this scenario, Silva 
et al. 14 suggest that the qx is a possible parameter that 


control the degree of anti-correlation between q from em¬ 
pirical rotational distribution and the cluster age t given 
by term q ss q»o(1 — Af/g^), where At = t — to is the 
difference between the logarithms of the cluster ages, and 


10 for a review). As mentioned by Epstein & Pinson¬ 
neault 161, first Reiners & Mothanty j9 argue that rota¬ 
tion is related to the magnetic field strength rather than 
the surface magnetic flux as adopted in Kawaler Jfi]. This 
dependence in R 2 modifies the structure of the angular 
momentum loss law, indicating that the Kawaler model 
presents a weak dependence on the stellar radius. And 
second, these authors assume that the saturation angu¬ 
lar velocity ui sa t does not depend on stellar mass. In this 


P-2 

















A nonextensive view of the stellar braking indices 


case, LOgat is insensitive to structural differences in stellar 
interior. In this context, Reiners & Mothanty [9] restrict 
the rate of angular momentum loss to radial fields and 
a dynamo law a = 1.5 for unsaturated regime, given by 
following equations 


d J 
d t 


= -C 


n 


(R 16 \ 1/3 

\AP ) 


( 1 ) 


furthermore, we assume that stars spin down as a solid 
body. 

Two pairs of equations were combined to obtain the 
time dependence of v. First, from eqs. (|3j) and the deriva¬ 
tive of ([5]), we have 


v(t) = Vo exp 


(t -1 0 ) 

Tl 


(to<t< t sat ) 


( 6 ) 


for f l > Lo sa t . Already, for f l < Lo sa t, we have 


— --r 
At 


= I (- 


1/3 


1/3' 


( 2 ) 


with C - 3 Kg2k , K1/ 

Here, G is the gravitational constant, ui sa t is the thresh¬ 
old angular velocity beyond which saturation occurs, B sat 
denotes the saturation field strength, Ky is a calibration 
constant and R, M and M represent the stellar radius, 
mass and mass-loss rate in units of 10 _14 MQyr _1 , respec¬ 
tively. 

For a future comparison with the parameters obtained 
in nonextensive Kawaler model presents in de Freitas & 
De Medeiros 10|, we eliminate the restrictions in a and 
in the magnetic field geometry N adopted by Reiners & 
Mothanty [9]. After a rather tedious algebraic manipula¬ 
tion, the generalized Reiners & Mothanty torque, for the 
different saturation regimes, where the parameters a and 
N are free, is given by 


dJ 

df 


= -C 


n 


r 1 


16—SN \ 1/(3—4JV)' 


M 2 


(3) 


for > ui sa t- In this regime, equation above is not depend 
on the parameter a. Already, for fI < ui aa t, we have 


dJ _ > 

At 


with C = 


Q \ 8a/(3—4iV) , Rl 
Usat) ^ V M 2 


16-8JV \ 1/(3—4JV)' 


(4) 


,4 JV 


C 3 


1/(3—4 V) 


If a = 1.5 and N = 0 

(later we will replace N = 2), we recovered the Reiners & 
Mohanty torque. 

In the present study, we assume that since the moment 
of inertia / and stellar radius R changes slowly during 
main-sequence evolution, the angular momentum loss law 
can be simplified by considering the limit of AJ/At for 
constants I and R, that is, this loss law is fully specified 
by the rotational velocity, while the star is braked by the 
stellar wind 17 . We also consider that in the absence 


of angular momentum loss, equatorial rotational velocity 
v can be determined by the simple condition of angular 
momentum conservation, denoted by 


j = ' L r 


(5) 


with 


C' ( r 16 ~ sn 

T ^ M 2 


1/(3—4AT)‘ 


— 1 


(7) 


for the saturated domain. 

On the other hand, by combining eqs. Q and the 
derivative of |5|, we obtain the equation 


v{t) = v sat 


1 tsat) 

T 


— (3—4JV)/8a 


it > t sat ) (8) 


with 


8a C' (R X3 ~ m 
3-ANT ^ M 2 


1/(3—41V)- 


— 1 


(9) 


In this respect, the t\ and r are the MS spin-down 
timescales in the saturated and unsaturated regimes, re¬ 
spectively. The other parameters, to, Vq and t sa t, are the 
age at which the star arrives on the MS, the rotational 
velocity at that time and subsequent age at which the 
star slows down to values below the critical velocity v sa t, 
respectively 9|. Time t sa t can be determined by setting 
v(t) = v S at in eq. ©>• Thus, we have 


t 


sat 


t 0 + Tl In 



( 10 ) 


Equation § reveals that in the saturated regime the 
star slows down exponentially from time t 0 up to velocity 
Vsat in time t sa t. Already, eq. ([8]) shows that from this 
time onward, the spin-down rate decreases to a power law 
given by —(3 — AN)/8a. Thus, the star remains within a 
factor of a few v sa t for the remainder of its MS lifetime. 


Nonextensive approach for stellar rotation-age 
relation for the unsaturated regime. Applying 
Tsallis statistics 18 , we can assume that eq. © follows 
a simple linear differential equation in the form 


d 

At 




from which the solution is given by 


( 11 ) 


v(t) = vq exp [-Ai (t - to)] • 


( 12 ) 


According to Lyra et al. [19 
nential behavior presented in eq. (12), eq. 


in contrast to the expo- 
must be 


replaced with a slower power law at critical points where 
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Fig. 1: Plot q from model against the field topology TV. Straight lines represent the Kawaler entropic index, while curves denote 
Reiners & Mothanty one. The vertical dashed dot line shows that exist only one valor for TV where q’s are equals, i.e., when 
TV = 2. The letter A represents the lowest value of qitM • Already, B denotes the Reiners & Mothanty [9] exponent, i.e., 

And C is the maximum possible value of q among the models. Plot also shows that the Skumanich law can not be represented 
on the generalized Reiners & Mothanty torque. 


long correlations develop. In this context, a similar expres¬ 
sion can be created to characterize the possible behavior 
of stellar rotation distribution in an unsaturated regime. 
As proposed by de Freitas & De Medeiros 10], the spin 


-f— 

d t \VsatJ 


= -A„ 


Vsat 


v(t) = V sat [1 + (q - 1) X q (t - t sat )} 1 


where [1 + (1 — g)#] 1 ^ 1 ^ = exp 9 (x) is known in nonex- 
tensive scenario as g-exponential (see 


scribed as a non-linear differential equation in form (13) 
with solution (14). When q = 1, eq. (12) is recovered, 


<ZitM can be described by a pair (a, TV) given by 

8 a 


<lRM = 1 + 


47V- 5’ 


(15) 


down law of a group of the stars limited by mass is usually 
given in the generalized form 


(\ > 0; g > 1), (13) 


if g = 3 the spin down is determined by torque generated 
by magnetically controlled stellar winds. 

Integrating the equation above, we have 


whereas for saturated regime, qRM = 1 due to the usual 
exponential function. As showed in table [l] the parame¬ 
ters a and TV are as a function of spectral type, i.e., the 
braking index q is related to stellar mass and evolution of 
stellar magnetic field. In Fig. [Tj we have the plot of index 
q against the field topology TV. Straight lines represent the 
Kawaler entropic index given by identity qx = 1 + 4a7V/3, 
while that the curves denote the Reiners & Mothanty pre- 


(14) 


scription given by term (15). For such model, figure shows 


20 for a review 


of g-calculus). In our context, X q indicates the braking 
strength. 

Thus, in the nonextensive scenario, eq. ® is well de¬ 


indicating that a system governed by a saturated regime 
is in thermodynamic equilibrium. By contrast, when q 
differs from unity, a system controlled by unsaturation is 
out of equilibrium. 

In agreement with our unsaturated regime model and 
replacing TV by TV — 2, we find that the “braking index” 


lines and curves as an iso-a (specifically a = 1, a = 1.5 and 
a = 2). The vertical dashed dot line shows that exist only 
one value for TV where qx and q^M are equals, i.e., TV = 2. 
From referred Figure, the letter A represents the possible 
lowest value of qRM- Already, B denotes the exponent 
adopted by Reiners & Mothanty i-e., qRM = 5. And 
C is the maximum possible value of q among the models. 
Figure also shows that the Skumanich law can not be rep¬ 
resented on the generalized Reiners & Mothanty torque. 

In general, the g-parameter is related to the degree of 
nonextensivity that emerges within thermodynamic for¬ 
malism proposed in Tsallis 11 . As revealed in eq. (15), 


qRM is a function of magnetic field topology TV and dy¬ 
namo law a, which depend on stellar evolution. According 
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Table 1: Best parameter values of our unsaturated model using 
eq. (141. The values of N using eq. (151 also are shown. 


Stars 

(M/M 0 ) 

IRM 

N a =2-,N a=1 

a N= o; a N=2 

F0-F5 

1.36 

1.80T0.1 

6.25; 3.75 

-0.5; 0.3 

F6-F9 

1.22 

1.96T0.2 

5.42; 3.33 

-0.6; 0.36 

G0-G5 

1.11 

2.18T0.2 

4.64; 2.95 

-0.74; 0.44 

G6-G9 

0.98 

4.38T0.6 

2.43; 1.84 

-2.11: 1.26 

All F 


1.91T0.1 

5.65; 3.45 

-0.57; 0.34 

All G 


2.18T0.4 

4.64; 2.95 

-0.74; 0.44 


to Kawaler 6| and de Freitas & De Medeiros 10], small 
N values result in a weak wind that acts on the MS phase 
of evolution, while winds with large N values remove sig¬ 
nificant amounts of angular momentum early in the Pre- 
Main-Sequence (PMS). In this phase, for a given mass, 
maximum rotational velocity v decreases as N increases. 
This result is significant because, within a thermostatisti- 
cal framework in which eq. ([8]) naturally emerges and for 
a given value of a and N, we obtain the scale laws found in 
the literature, such as those proposed by Skumanich |3j, 
Soderblom et al. 21 and Reiners & Mothanty 9 . 


CONCLUSIONS. — In summary, our study presents 
a new approach on the angular momentum loss law found 
in literature. The present analysis shows that the gen¬ 
eralized torques of Kawaler and Reiners & Mothanty are 
clearly antagonistic, except for N = 2. Table [l] demon¬ 
strates clearly that the modified Reiners & Mothanty 
model is valid only for spectral types beyond G6. Accord¬ 
ing to Fig. [l] the lowest valor of qrm is 3.67. As a result 
from this figure, the Skumanich law can not be represented 
on the generalized Reiners & Mothanty relationship, i.e., 
this model does not explain the rotational behavior of solar 
type stars. This limitation is mainly due to sensitivity on 
the stellar radius given by term I? 16 / 3 , i.e., the angular mo¬ 
mentum loss rate would have become entirely dependent 
on the stellar radius. In Kawaler model, this dependence 
is given by maximum term R 2 , where depends entirely 
on the stellar velocity for unsaturated stars. However, our 
results reveal which Reiners & Mothanty model is in good 
agreement for low-mass stars, specifically, M-type stars. 
In other words, the Reiners & Mothanty model is a good 
approximation for fully convective stars. 

As mentioned by de Freitas & De Medeiros 


10 , the en- 


tropic index q is an indicator of the memory of system or, 
in our case, the memory of primordial angular momentum. 
As q increases from F to M, very low-mass stars maintain 
this type of memory. In fact, the memory grows at same 
scale of q. In contrast, because the much shorter interac¬ 
tion of the young star with its circumstellar disc, higher 
mass stars do not maintain the memory of initial angular 
momentum for long time. According to de Freitas & De 
Medeiros 10 and our present results, both nonextensive 


models show that the Kawaler framework is compatible 


with a wider mass range, while the Reiners & Mothanty 
model is restricted to low-mass stars less than G6 masses. 
This result suggests that the decrease in mass would lead 
to weaker power law than the Skumanich law. Finally, the 
results point out that a strong dependence of angular mo¬ 
mentum loss rate on stellar radius is limited rather to very 
low-mass stars encountering discrepancies when applied to 
F and G stars. 
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